Socioeconomic Networks with Long-Range Interactions 
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We study a modified version of a model previously proposed by Jackson and Wolinsky to account 
for communicating information and allocating goods in socioeconomic networks. In the model, 
the utility function of each node is given by a weighted sum of contributions from all accessible 
nodes. The weights, parametrized by the variable 5, decrease with distance. We introduce a growth 
mechanism where new nodes attach to the existing network preferentially by utility. By increasing 
5, the network structure evolves from a power-law to an exponential degree distribution, passing 
through a regime characterised by shorter average path length, lower degree assortativity and higher 
central point dominance. In the second part of the paper we compare different network structures 
in terms of the average utility received by each node. We show that power-law networks provide 
higher average utility than Poisson random networks. This provides a possible justification for the 
ubiquitousness of scale-free networks in the real world. 



I. INTRODUCTION 

The study of socioeconomic networks is a burgeoning 
field in the physics and economics literature, with ma- 
jor progress having been attained over the last decade 
HI) Hi H) 0) H) @| • Individuals and firms interact through 
networks to share information and resources, exchange 
goods and credit, make new friendships or partnerships 
etc. The structure of the network through which interac- 
tions take place may thus have an important effect on the 
success of the individual or the productivity of the firm 
Furthermore, the network of interactions among so- 
cioeconomic agents plays an important role for the stabil- 
ity and efficiency of socioeconomic systems Q- Theories 
about how interaction networks form are thus essential 
for a deeper understanding of the development and orga- 
nization of society as a whole. 

The economics literature focuses mainly on equilib- 
rium networks and the network formation mechanisms 
are based on utility maximization and costs minimiza- 
tion. The aim of most economic papers is to identify, 
among the set of equilibrium networks, the geometry that 
optimizes efficiency [36| in the sense of social benefit. 
Likewise, economists are interested in the stability |37j 
of equilibrium networks under link deletion, addition or 
rewiring [J, . A shortcoming of these models is that the 
equilibrium networks are often too simple in their geom- 
etry (stars, complete networks, interlinked stars, etc.), 
typically as a consequence of the symmetries that need 
to be assumed in the payoff functions in order to make 
the models analytically tractable Q • 

The physics literature, instead, has mainly focused on 



the characterization of the structure of real networks and 
proposed dynamic models, mostly based on probabilistic 
rules, capable of reproducing the observed geometrical 
structures (Poisson, stretched exponential and scale-free 
networks) gU^,!!!- 

In this paper we try to combine the physics and 
economic approaches, by introducing a stochastic net- 
work formation mechanism inspired by economists' util- 
ity maximization models, which naturally extends the 
well known physicists' preferential attachment rule [l^ . 

One of the most interesting models of socioeconomic 
network formation was introduced by Jackson and Wolin- 
sky in 1996 [l]. In their model, the formation and evolu- 
tion of links is driven by a utility maximization mecha- 
nism. The model is based on the assumption that agents 
may derive benefit not only from the nodes to which 
they are directly connected (their nearest neighbours), 
but also from the ones they are connected to indirectly 
(possibly via long paths). Less distant connections are 
more valuable than more distant ones, but connections 
to the nearest neighbours are costly. The utility of node 
i is defined as: 



(1) 
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where the contribution to the utility of i from j may 
depend on the weight Wij of the edge between i and j 
(or, alternatively, on the fitness of node j); < 5 < 
1 captures the idea that the utility gain from indirect 
connections decreases with distance; dij is the number 
of links in the shortest path between i and j {dij — oo 
if there is no path between i and j); V{i) is the set of 
nearest neighbours of i; and Cy are the (node specific) 
costs to establish a directed connection between i and j 
[3^ . Costs can also be differentiated in costs of initially 
creating or maintaining an edge 
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The papers by Jackson and Wolinsky [l| , as well as the 
one by Bala and Goyal [3], are mainly concerned with 
stability and efficiency of the network resulting from dif- 
ferent dynamic updating rules. In particular, Jackson 
and Wolinsky study pairwise stability when agents can 
only update a link at a time (either delete it or create 
it), while Bala and Goyal allow agents to rearrange all 
their connections at once. The updating is determinis- 
tic in both models, and a new configuration is accepted 
only if it increases the utility of the agent. These two 
papers show that the star network is both efficient and 
stable for a wide range of the parameters when S — 1. 
Nonetheless, a multiplicity of network architectures exist 
in [2| for < (S < 1 which could be a strict Nash equilib- 
ria, and to which the system may converge depending of 
the initial conditions. Feri 13] has shown that for suffi- 
ciently large networks the star network is stochastically 
stable for almost all the range of parameters, even for 
0<6<1. 

Here we focus on the connections model of Jackson and 
Wolinsky, i.e. the case Wij = 1, Wu = and c.y = c. In 
this case, the utility can be written as 

u^Jf: E E c=E5'zP-cz« (2) 

1=1 {k\dik=i} jevii) 1=1 

where the sum in / is over all shortest paths of length I 
from node i, the sum in k is over all nodes whose shortest 
path from i is dik = /, Imax is the path length of the node 
the furthest away from node i, and is the number of 
Zth- nearest neighbours of node i. The utility of a node 
is expressed in ([2]) as a weighted sum of the number of 
nodes accessible from i on outward " layers" of increasing 
distance from i. Thus, we start at node i and multiply 
S by the number of nodes that are joined by an edge to 
i-this being the first layer. We then add S"^ times the 
number of nodes that are joined by an edge to a node 
in the first layer-this is the second layer. We continue in 
this way until no new nodes are found. Hence, expression 
([2]) incorporates implicitly the well known breath-first 
search algorithm ^14| . In this paper, we will consider 
only networks with zero costs. Therefore, equation ^ 
becomes: 

= E -^'4^^ (3) 
1=1 

In the first part of the paper we focus on a specific 
network growth mechanism and examine the resulting 
network topology. If each new node attached determin- 
istically to the existing node with maximal utility, the 
resulting network would be a star. The randomness gen- 
erated by a probabilistic attaching rule can be interpreted 
as costs and barriers to gather information, or bounded 
rationality, all of which limit the ability to establish links 
in an optimal way, thus possibly generating more realis- 
tic geometries than the star network. It is thus, worth- 
while to ask which network topologies are to be found 




FIG. 1: Schematic layout of network growth when m = mo = 
1. The addition of a new node, 6, implies an increase of the 
utility of nodes 1 to 5 which is simply S"^, where d is the 
path length from node 6. The simplicity of this updating 
mechanism allowed simulations to be run with N = 10^ when 
m = 1. 



when new nodes arrive steadily and create links with ex- 
isting nodes in a probabilistic way, proportionally to the 
utility of existing nodes. In this way, we build on the 
preferential attachment growth rule of Barabasi and Al- 
bert 0, [l2| which can be recovered from equation ^ 

(i) 

when Imax = 1- Furthermore, preferential attachment is, 
arguably, the most extensively studied mechanism of net- 
work formation and one that has revealed insights into 
properties observed in real networks. Therefore, it is im- 
portant to understand the robustness of the specific rule 
of linear preferential attachment by node degree, which 
is one of the aims of this paper. 

Often, the specific network growth mechanisms are un- 
known and only the topology of the equilibrium network 
can be extracted from data. One obvious question is then 
how the observed equilibrium networks rank in terms of 
their efficiency, e.g. Erdos-Renyi versus scale-free net- 
works. We address this question in the second part of 
the paper and derive analytical results by using the gen- 
erating function approach |,15i] . We show that power-law 
networks are more efficient than Poisson random network 
when individual utility is defined by ([3]) thus providing 
a possible explanation for why scale-free networks are so 
ubiquitous. 



II. GROWING NETWORKS 

In the classic Barabasi and Albert model [13], a net- 
work is grown by adding, at every time step, a new node 
that attaches to m existing nodes with a probability pro- 
portional to their degree, Il{ki) = ki/^^^-^ kj. At time 
t, the resulting network has size Nt = toq -I- t, where 
Too > to is the size of the (fully connected) network at 
time t — 0. Preferential attachment models were in fact 
introduced in the literature already by Yule and Si- 
mon [H]. 

The preferential attachment mechanism generates a 
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FIG. 2: (Color online) Cumulative distribution function of degree (panels a), b) and c)), utility (panels d), e) and f)) and 
betweenness centrality (panels g), h) and i)) for several values of 5 £ ]0, 1], and m = 1, 2 and 5. We also plot the corresponding 
distribution of degree and betweenness for the BA model (curves were shifted vertically). Simulations were averaged over 30 
runs in networks with A*' — 10^ (m = 1) or A = 5 x 10^ (m = 2 and 5). Coloured bands around the curves are 95% confidence 
intervals. 



scale-free probability density of incoming links that leads 
to the stationary result p{k) — 2m? /k^ , with 7 = 3 in- 
dependently of m. The model is also characterized by 
a clustering coefficient larger than the one found for the 
Erdos Renyi networks (for m > 1) and no clear assorta- 
tive/disassortative behaviour [9|]. 

Several models have been proposed lately to investi- 
gate extensions of the preferential attachment mechanism 
through edge removal and rewiring, inheritance, redirec- 
tion or copying; node competition, aging and capacity 
constraints; and accelerated growth of networks to name 
just a few (see H,iiQ., 17, 35] for reviews). 

Of particular relevance to our approach are fitness 
models [IMISIISIIH) where the probability of attaching 
to a node is proportional to node fitness 



fiki 



(4) 



Here we extend the preferential attachment rule by in- 



troducing a growing mechanism inspired on the work of 
Jackson and Wolinsky Our contribution is to pro- 
pose preferential attachment by node utility. Thus, the 
probability that a new node j will be connected to an 
existing node i depends on the utility of i, such that 



n, 



(5) 



where the utility of node z, Ui, is given by Attach- 
ment happens with uniform distribution for (5 = 1 gen- 
erating an exponential distribution of node degree [39|. 
While the model is not defined for (5 = 0, the preferential 
attachment rule (O is invariant up to multiplicative fac- 
tors in so for (5 7^ the qualitative behaviour of the 
model remains unchanged if we define utility as 



-1 Ji) 



(6) 



1=2 
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FIG. 3: Kamada-Kawai spring layout ^16^ for m = 1 and — 10^. Panel a) is a sample layout for S — 0.01, b) for 5 = 0.2, 
c) for S = 0.7 and d) for 5 = 1. On each panel, nodes are coloured by their utility on a gray scale from minimal (white) to 
maximal (black) utility. 



where ki is the degree of node i. Thus, as (5 — > our 
model converges to the Barabasi-Albert model and gen- 
erates a scale-free network. 

Our model has resemblances with the fitness models 



discussed above. However, there is a fundamental dis- 
crepancy: we regard utility as a time-dependent measure 
of node fitness, whereas existing models assume that node 
fitness does not change with time. 
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FIG. 4: (Color online) Plots of a) average path length, b) central point dominance, c) assortativity coefficient, and d) normalized 
maximum degree as a function of 5 for the simulation results when m = 1, 2 and 5. Curves in panel a) were scaled by the values 
at (5 = 1 (exponential network). The BA values are indicated by the corresponding full symbols. Simulations were averaged 
over 30 runs and coloured bands around the curves are 95% confidence intervals. 



At each time step, a new node j joins the network and 
the utility of existing nodes changes. When m = 1, the 
utility increment to an existing node i at distance / from j 
is given by Aui — S' and therefore, at each time step, the 
computation of utihty for the network can be completed 
in 0{N) time. Figure[l]is a diagram of a possible network 
configuration with m — toq = 1 after t = 5 time steps, 
showing the change in utility of existing nodes, Aui = 
SK When m > 1, the increment in the utility of node 
i depends on the existing network geometry and Aui > 
5^. Therefore, when m > 1, we need to re-compute the 
utility of all existing nodes at every time step, and the 
computation runs in 0{N'^) time as it involves running 
a breadth-first-search algorithm from every node. This 
is the reason why we have ran simulations for N — IQ^ 
when m = 1, but only up to A'^ = 5 x 10'^ when m > 1. 

Existing nodes i at a higher distance than a certain 
^max from new node j receive a contribution Aui = 



Sd{j,i) ^ -^Q-precison iggg ^j^^n the number of 

significant digits that the computer can store (typically 
precison = 32 in double precision), and do not need to 
have their utility updated in the simulations. This max- 



> 



precision 
logio S 



(7) 



Our implementation of the algorithm updates the utility 
of all nodes accessible from the new node j up to distance 
^max = —32/ logj^o ^- The code was implemented in C++ 
and ran on a Condor framework (high throughput com- 
puting) [i^l for several values of S. Ensemble averages 
were taken over 30 runs. 

Expressions ([3]) and ([6]) predict the existence of two 
distinct regimes: a scale-free regime as (5 ^ (<5 7^ 0) and 
a random growth regime for which degree distribution is 
exponential at i5 = 1. We are interested in exploring how 
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FIG. 5: (Color online) Plots of average k/kmax and u/umax 
as a function of the entry time for the first 100 nodes when 
A'' = 10"" and ra = 1. Simulations were averaged over 30 runs. 



the network evolves from one limit regime to the other 
as we increase 5. 

Figure d^-c) shows the distribution of degree for m — 
1,2 and 5. We also plot the corresponding distribution 
for the BA model (solid curves shifted vertically). For 
very small 5 ~ 0.01, preferential attachment by degree is 
indistinguishable from preferential attachment by utility 
and the probability distribution of both quantities decay 
like p {x) ~ with 7 = 3. The power-law decay in the 
BA model is known to be a peculiarity of the linear pref- 
erential attachment mechanism and is destroyed by small 
perturbations like, for example, a non-linear attachment 
rule n(fci) ~ kf [13] ■ Here we also observe a depart from 
the scale-free regime as we increase 5. Furthermore, in 
the Barabasi-Albert model the degree distribution de- 
cays as a power-law with exponent 7 = 8 independently 
of m. In our model, increasing m has the effect of ho- 
mogenizing the utility of the nodes (the distance between 
pairs of nodes decreases with increasing connectivity in 
the networks). Consequently, deviations from the power- 
law decay are observed at lower and lower values of 5 as 
we increase m. 

Betweenness centrality is plot in Figure as m is 
varied. Recent results have shown that the distribution 
of loads (or betweenness) scales with a power law [2^, [23] 
p{g) ^ where a = 2 for a tree (and hence for m = 1). 
This justifies the collapse of the curves of the distribution 
of betweenness in Figure [5^). As can be observed in Fig- 
ure [2h) and i), the distribution of betweenness deviates 
from the power-law behaviour as m is increased. 

For intermediate values of 6 and m = 1 a number of 
interesting features appear. First of all, we observe that 
the utility distribution becomes step-like (Figure[2Ji) sug- 
gesting the presence of subsets of nodes that share similar 



utilities. 

This phenomenon can also be inferred from the net- 
work layouts in Figure [3] (for networks of 10^ nodes with 
m = 1), which are produced using the Kamada-Kawai 
spring layout [16]. Essentially, the Kamada-Kawai lay- 
out assigns stronger springs to vertices that are closer in 
the graph-theoretic sense (i.e., by following edges) and 
therefore places them closer together. In the case m = 1 
(a tree), nodes close to the hubs on the graph layout 
will also be close in graph terms and therefore we can 
interpret the layout heuristically: these nodes have sim- 
ilar utility. When 5 = 0.01 (close to the BA scale-free 
regime), the layout shows a few utility hubs (the dark 
vertices) surrounded by clouds of nodes that disperse as 
we move further away from the hubs; for ^ = 0.2 denser 
clouds of nodes cluster around a smaller number of hubs, 
and can still be observed farther away from the hub; for 
higher 5 the clouds start dispersing; eventually for 5 = 1 
all nodes have the same utility. 
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FIG. 6: (Color online) Average utility for the simulation re- 
sults (solid curve and symbols) and the BA model (dashed 
curve, open symbols) for m = 1 (21 = 2), 2 {z\ = 4) and 5 
(z\ = 10) when N = 5 x 10^. We also plot the average utility 
of a star with N = 5 x 10^ nodes. Simulation results were 
averaged over 30 runs. Coloured bands around the curves are 
95% confidence intervals. 

The rearrangement of the network as S increases from 
zero gives rise to non monotonous behaviour of the av- 
erage path length I, degree assortativity r degree , and 
central point dominance [2H ] . Both I and r degree show 
a minimum for S 0.15 {N = 10^ and m = 1) and C'g 
has a maximum around the same value. This behaviour 
can be observed in FigureH^-c). The average path length 
I((5) is measured relative to the path length of a random 
growing network (i.e. relative to 1(1)). The scale-free 
regime is characterised by a shorter path length than the 
random growth regime. Here we observe an even further 
contraction of the network for values of ^ up to ~ 0.5. 
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Note that the average path length of a star network is 
I* = 2 for N large. When normahzing with — 20.16 
we get a value of I* ^ 0.1. While this value is still much 
smaller than 7(0.15), the network contraction seems to 
indicate a move toward a more star-like configuration. 
To further investigate this point, we compute the the 
central point dominance measure introduced by Freeman 
[26] and plot it in Figure [Us) . This measure is defined 
as the average difference in node centrality (measured by 
node total betweenness) between the most central point 
and all the others. The central point dominance takes a 
value between zero (for a graph where all points have the 
same centrality) and one (for the wheel or star graph). 
The maximum of C'^ for J ^ 0.15 in Figure IUd) confirms 
that the network is becoming more star like around these 
values of S. 

Next, we plot the assortativity of the network in Figure 
SI:) . We implement as measure of assortativity the degree 
assortativity r degree [25| which takes values from —1 to 
+1: negative values for disassortative networks, if the 
networks are neither dissasortative or assortative, and 
+1 for fully assortative networks. This value approaches 
zero for large N in the BA model [13] and is negative for 
a star. Our model shows a lower assortative mixing than 
the BA model for values of 5 up to ~ 0.5. The decrease 
in Tdegree IS also consistcnt with our hypothesis that the 
network is becoming more star like at 5 ~ 0.15. While 
the network goes through this rearrangement the degree 
of the most connected node (as a fraction of the total 
number of links) is nonetheless monotonously decreasing 
with S, as shown in Figure UJi). The same behaviour is 
observed for the utility of the most connected node (not 
shown here). This reveals that as new nodes are added 
to the network, they do attach on average closer to the 
hubs as S increases, generating a more compact network, 
but not directly to them. 

To gain further insights into the structural changes 
that take place as S increases from zero, we analyze the 
role of entry time on node connectivity by computing 
the ratios u/umax and k/kmax for the first 100 nodes as 
a function of the time at which they entered the network. 
The ratios plotted in Figure[S]are averaged over 30 differ- 
ent simulations, for networks with N = 10^ and m — 1. 
The plot shows that for any 5 the initial node is likely 
to acquire the highest fraction of links and utility. More- 
over, for 5 = 0.15 both the degree and utility ratios decay 
faster than in the scale-free regime. This reveals that, for 
5 ^ 0.15, the earlier nodes receive both a higher relative 
degree and utility than in the scale-free case. In other 
words, the earlier nodes are stronger hubs than in the 
scale-free case and thus the network arranges in a more 
star like configuration around i5 = 0.15. As 5 increases 
further, the slope of the utility ratio becomes lower than 
in the scale-free regime. In this range, the utility differ- 
ences between old and new nodes are not large enough 
to create well defined utility or degree hubs. 

Figure [5] also highlights the fundamental mechanism 
of structure formation in our model: the fine relation be- 



tween degree and utility as i5 is varied. In the scale-free 
regime, preferential attachment by utility is equivalent to 
preferential attachment by degree and node degree and 
utility assume the same value as the network grows. At 
6 ^ 0.15, we observe a gap between the two scaled quanti- 
ties, suggesting a discrepancy between degree and utility 
for the higher order neighbours of the utility hubs. This 
gap is larger for 5 = 0.15 than 5 = 0.4, indicating that 
around 5 ~ 0.15 the growth mechanism is generated by a 
variable (node utility) which is considerably independent 
of node degree and thus revealing why this is the region 
where the network displays more interesting structure. 
As 5 increases towards 1, the influence of random net- 
work growth becomes more important and this structure 
generating mechanism disappears. 

Finally, we investigate how average node utility com- 
pares in networks generated with our preferential utility 
attachment, the scale- free regime (here generated via the 
BA preferential attachment mechanism) and a star net- 
work. The average utility of a star network is given by: 

il*(<5) = fci (^l + <5^^^ (8) 

where zi = 2{N — 1)/N. For N large, zi ~ 2 and 
M* (d) ^ NS'^. In Figure [6] we plot the average utility for 
networks in our model at different S when = 5 x 10'^ 
and zi = 2,4 and 10 and compare that to the correspond- 
ing average utility in the BA model and a star with the 
same N. The plot shows that the BA scale-free network 
has a higher utility than the network generated via our 
preferential utility mechanism at all values of delta and 
zi. Networks in our model become more star-like for 
6 0.5, but this implies an increase in the utility of only 
a small number of nodes (the early nodes). Therefore, 
the average utility of nodes in the network is still higher 
in the BA model than in our model for the same values 
of S due to the scale-free structure of the former. Com- 
parisons with the star network can only be made when 
^1 = 2 as this is the average degree of the star network 
when N is large. Figure [6] confirms that the star network 
has the highest utility for this value of zi among all the 
networks we study. Nevertheless, the star network can 
only be achieved if agents are perfectly rational and have 
access to full information (in which case the attachment 
mechanism would be deterministic). This is rarely the 
case in real word situation, thus the comparison with the 
star is of little practical relevance. 

III. ANALYTICAL RESULTS FOR RANDOM 
NETWORKS 

An interesting question to ask (for example, from the 
point of view of the social planner) is how network topolo- 
gies rank against each other and which network struc- 
ture maximizes the total, or the average utility (networks 
that satisfy this condition are said to be efficient in eco- 
nomics). We show that it is possible to derive analytical 
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FIG. 7: (Color online) Average number of first and second- 
neighbours {zi (7, a) and 22 (7, a)) in networks with degre dis- 
tribution given by (If 5|) . From left to right, we plot zi (7, a) 
(full curves) and Z2 (7, a) (dashed curves) for a — (black), 
f (blue), 2 (green) and 3 (red). The values of zi (7 < 3, a) 
for which Z2 (7, a) is not defined are plot in grey, as well as 
regions of the curves for which zi (7,0) > 22 (7, «)• The cir- 
cles denote the intersection of the two curves, zi (7, a) and 
Z2 (7, a), for each value of a. 
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FIG. 8: (Color online) Analytical results for average utility 
in networks with power-law (full curves) and Poisson (dashed 
curves) degree distributions as a function of J, z\ and 7 for 
= 10^. Curves have been shifted vertically for different 
values of 7 for clarity. Values of z\ increase from bottom to 
top. 



results for the average utility in Poisson and power-law 
networks. By comparing average utility in different net- 
work topologies with the same size and the same average 
degree, we show that power-law networks are more ef- 
ficient than Poisson random networks (even though less 
efficient than the star) . The effect of costs on u (5) is a 
constant term for all networks with the same N and z\. 
Therefore, without loss of generality, we choose c = in 
the analysis below. 

If the sum in ([3]) was to be evaluated up to distance 
l^max = 1 for every node, expression ([9]) would simplify to 
u{S) = 5zi^ i.e. average utility would be independent of 
the specific network topology and all networks with the 
same number of nodes and links would be equally effi- 
cient. Thus we need to introduce long range interactions 
{Imax > 1) to be able to rank networks in terms of their 
efficiency. 

To derive an expression for average utility in generic 
random networks with N large, we average both sides of 
®: 



zi 



(9) 



1=1 



where zi is the avera ge n umber of ^th neighbours of a 
node. Newman et al. [28l [2^ define I via the expression 



closed form of expression ^ if we have access to analyti- 
cal expressions for I and z; . This can be accomplished by 
generating functions, which arc particularly useful when 
determining means, standard deviations and moments of 
distributions [isj . 

The average number of neighbours (average degree) 
and the average number of second neighbours of a node 
can be derived from the probability generating function 
of node degree. Go (x) = Y^^^oPkX^ ^ long as the de- 
gree distribution, p^, is specified. The beauty of the gen- 
erating function formalism is that one can derive zi as a 
function of zi and Z2 only [2^ [2^ [sO] : 



zi 



n '-1 



(11) 



Replacing equation (fTT|) in equation ^ yields 



X X^ix^V-^ ^zi{{SZy - 1) 
■.(d) = dzi2_^{dZ) = — — (12) 



1=1 



SZ - 1 



where Z = zijzx. For Z > \ and > -I- 1, which are 
conditions satisfied by most networks, I can be calculated 
as a function of A^, z^ and 22 from PH)) and pT|) as [28j : 



1 + ^ zz = Af 



(10) 



Now that we have expressed average utility in terms 
of the breadth-first search algorithm, we can derive a 



ln[(Ar-l)(Z-l)/zi + l] 
ln(Z) 



(13) 



In what follows, we investigate the behaviour of (I12p for 
Poisson and power-law random networks. 



9 



Poisson random networks are characterized by zi = 
pN and Z2 = [28], thus ^ yields = 

In ( ''^^^ ) / In(zi). In this case, becomes: 

up{N,S,z^) = ^ ^ (14) 

for N > zi + 1, < 6 <1 and zi > 1. 

Next, we consider power-law networks with degree dis- 
tribution of the form: 



Pk (7, a) 



1 



-{a + ky'^ ,a>0 (15) 



C(7,l + a) 

where the normalizing factor ^(7, a+1) = 



k) ^ is the Hurwitz zeta function (7 > 1). The gen- 
erating function for the probability distribution is given 
by 



Go(.,7,«) = EP^-^-^?fe^^ (16) 



fc=i 



C(7,a+1) 



where $(x,7, a) = ^'kLo (a+k)-' Lerch transcen- 

dent. For our purposes, only the first two derivatives of 
$(a;,7,a -I- 1) with respect to x are relevant as the av- 
erage number of first and second- neighbours are given, 

respectively, by zi (7, a) = 
. Hence 



and Z2 (7, a) 



d''Go(x) 



zi 7,a = — ,7>2Aa>0 17) 

C(7,a+ 1) 

Z2 (7, a) = ^^77 — -^-TV^^i ("^ - 1, a) - (a + l)zi (7, a) , 7 > 3 A a > 



(18) 



Thus 



C(7-l,a+l)^i (7-1, a) 



C(7,a-fl) zi(7,a) 
Substituting ^ and (HH) into ([121), we find 



-a-1 ,7>3Aa>0 



(19) 



In 



(a+2)(jV-l) zi(7-l,a)C(7-l,a+l)(jV-l) , l 
2l(7.a) 2i(7,a)^C(7,a+l) ^ 

V zi(7,<i)C(7.a+l) 



, iV > zi (7, a)-t-lA7 > 3Aa > OAZ (7, a) > 1 (20) 



and thus average utility is given by 

6zi (7, a) (^{5Z (7, a)) 



lsF(N,^,a) 



usF {N, S, 7, a) = 



6Zi^,a)-l 



,iV>zi(7,a) + lA0<(5<lA7>3Aa>0AZ(7,a)>l 



where 21(7,0), Z{j,a) and lsF{N,^,a) are given by 
(fT7| . (HH) and respectively. 

When a = 0, the distribution of degree, (fT5|) . becomes 
a pure power-law pfc (7) = ■^^k^''. In this case, we have 

C(7,a+l)la=o = C(7) and 'i>(.T, 7, a + l)L.o = 
therefore (fTBl) becomes 



^i(7)L=o-^^^ ,7>2 (23) 

C(7-2)-C(7-l) .o,^ 
^2 (7)L=o = ' 7 > 3 (24) 

Therefore, in pure power-law networks, when N —> 00, 
the average number of second- neighbours, 22 (7)7 is fi- 
nite only for 7 > 3. However, the Riemann zeta func- 
tion, C (7), is a decreasing function of 7 (for 7 > 3) and 
zi (7 = 3) = 7rV6C(3) ~ 1.36843. In other words, the 
This generating function is also obtained for the power- existence of za (7) implies zi (7) < zi (7 = 3) =i 1.36843, 
law distribution with exponential cut-off, proposed in (Sll, which is a non-realistically low value for average degree 
|32|, Pk (7, k) = Ck-^e~^/'', in the limit k ^ 00. 
Expression ([^^ implies 



Go (a;, 7) 



Li^ (a;) 

T(7r 



(22) 



in real networks. This explains why we have chosen the 
modified power-law distribution psp . 
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The generating function ITBl) encapsulates all the mo- 
ments of the degree distribution 28]. Hence, the ex- 
pressions for zi{-f,a) and 22(7,0), (fTT]) and p8|l . are 
only exact in the limit N — > 00. Further, Isf iN,j,a) 
and ugi? (iV, (5, 7,a), both of which depend on Z2 (7,0), 
are only defined where Z2 (7,0) is finite, i.e. for 7 > 3. 
Therefore, it is essential to understand the behaviour of 
zi (7, a) and Z2 (7, a) in power-law networks. Figure [7] 
shows zi (full curves) and Z2 (dashed curves) within the 
range 7 > 3 A Z > 1 (where Isf {N, 7, a) is defined) for, 
from left to right, a = 0, 1, 2 and 3. 

Having deduced closed-form expressions for average 
utility in Poisson and power-law networks, we can now 
compare both networks under the condition that zi is 
the same. Figure [8] is a plot of average utility versus 6 
when zi = {2, 4, 10} and N — 10^ for Poisson and power- 
law networks. The average utility of Poisson networks is 
completely specified by N, d and zi , but power-law net- 
works defined by (jlSp have one extra degree of freedom 
in zi (7,0). In this case, we compute zi numerically by 
solving (Hil) for zi (7, o) = {2, 4, 10} when 7 = {3.1, 4, 5}. 
For all cases studied power-law networks are more effi- 
cient than Poisson networks. 



IV. DISCUSSION 

The growth mechanism we have proposed is a natural 
extension of the Barabasi- Albert preferential attachment 
by degree to preferential attachment by node utility. Our 
analysis shows that for small values of S, the utility decay 
parameter, the network retains a scale- free structure that 
is nonetheless destroyed when S increases. We have iden- 
tified a regime in 6 where the network is characterized by 
a lower average path length and assortativity coefficient 
and a higher central point dominance than the scale-free 
network. In this regime, the distribution of utility is a 
step-like function and the network has a more star-like 
structure. 

The derivation of analytical expressions for average 
utility in Poisson and power-law networks reveals that 
the latter have higher u for the range of parameters 
that is of significance in real- world networks (zi > 2). 
This suggests a novel mechanism which may explain the 
ubiquitous presence of power-law networks, in partic- 
ular in situations where collaboration, interaction and 
information sharing among the nodes are of paramount 
relevance. 



Social networks are highly volatile. Friendships can 
be stable for a long time but occasional encounters may 
lead to creation of links that are never used again in the 
future. A dynamical model of network formation, where 
links not only can be heterogenous, but whose weights 
can change continuously over time, would be a more ap- 
propriate way to describe social interactions. Our as- 
sumption that Wij = 1 for all links is obviously a first or- 
der approximation and preferential growth, with no link 
rearrangements, is a crude description of social network 
formation. Nonetheless even this simple mechanism can 
highlight surprising features of the models (like in our 
case a smaller network diameter for intermediate values 
of S) and as such it is worth to investigate in more general 
contests than the original BA model. 

We have also assumed that the connection costs in our 
model are zero. This assumption is justified by the fact 
that if costs are node independent they do not play any 
role in the growing model. Similarly, costs do not play a 
significative role if we restrict the comparison of average 
utility in section IIIII to networks with the same size and 
the same average degree. Nonetheless in a more realistic 
model, where links can be rearranged over time, costs 
would also play an important role in determining the 
shape of the network. 

Further analysis taking into account both these effects 
is currently under development. 
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